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Abstract. New integrable lattice systems are introduced, their different integrable dis- 
cretization are obtained. Backlund transformations between these new systems and the rela- 
tivistic Toda lattice (in the both continuous and discrete time formulations) are established. 



1 Introduction 



This paper is devoted to integrable equations of the classical mechanics. More precisely, 
we shall deal here with equations of motion in the Newtonian form. 

We introduce here two new integrable continuous time lattice systems, and present 
several novel integrable discrete time systems. 

The first continuous time system is: 

x k = x k+1 exp(x fe+ i - x k ) - exp(2(x A;+ i - x k )) - exp(x k - x k -i) + exp(2(x fc - x fc _i)). 

(1.1) 

The second one: 

'x k = x 2 k (x k+1 exp(x k+1 - x k ) - Xfc_iexp(x fe - (1.2) 

To the author's knowledge, these systems have not appeared in the literature, despite 
their beauty and possible physical applications. However, they are in a close relation to 
another well known integrable lattice - a relativistic Toda lattice: 

.. _ . . g 2 exp(x k+1 - x k ) . . g 2 exp(x k - x k _i) 

%k — X k+ iX k - — ■ r — X k X k ^\- — ■ r. 

l + g 2 exp(x fc+ i -x k ) 1 + g 1 exp(x fc - x k ^ ) 

More precisely, we shall show that fll.lf) and (|1.3| ) are connected by a sort of Backlund 
transformation, and the same is true for ( |1.2| ) and ( |1.3| ). 

We now write down integrable discretizations we propose for the lattices fll.lj) , ( |1 . 2|) . 
In the difference equations below x k = x k (t) are supposed to be functions of the discrete 
time t £ hZ, and x k = x k (t + h), x k = x k (t — h). Discretization of the first lattice (|1.1|): 

exp(x fc - x k ) - exp(x k - x k ) = 



\-hexp(x k+1 -x k ) + - --hex-p(x k -x k _i). (1.4) 



1 - hexp(x k+i -x k ) " """ 1 - hexp(x k - x k ^ x ) 



Discretization of the second lattice ( |1.2|) :: 
h h 



exp(x k - x k ) - 1 exp(x k - x k ) - 1 

exp(x fc+ i - x k ) - exp(x fc+ i -x k ) - exp(x fe - x k ^) + exp(x fc - £ fc _i). (1.5) 



1 



The same Backhand transformation as for the continuous time systems relates these 
systems of difference equations to the discrete time relativistic Toda lattice: 



exp(x k - x k/ 



1 (l + g 2 exp(x k+1 - x k )j (l + g 2 exp(x k - x k -i) 



exp(x k - x k ) - 1 (l + gi exp{x k+1 -x k )) (l + g 2 exp{x k - x k _ ± ) 



1.6) 



A modification of the construction leading to the discrete time systems above allows 
to derive several further nice discretizations. For example, for the lattice ( p..l|) : 

(l + hexp(x k+1 - x fc ))(l - hexp(x k+1 -x k )) 

exp{x k - 2x k + x k ) = T -y-r z -y, (1.7) 

(1 + hexp(x k - x k -i)j[l - hexp(x k - x k -i) 



and for the relativistic Toda lattice ( |1.3|) : 

exp(-x k + x k ) - 1 _ i 1 + ^ 2 exp(x fc+ i - x k yj (l + g 2 exp(x k - x k ^ 



exp(-x k + x k )-l (l + g 2 exp(x k+1 -x k )) (l + g 2 exp(x k - x k ^) 



'U 



(The last system resembles very much the previous discretization of the relativistic Toda 
lattice ( |1.6| ), however the relation between them is far from trivial). 

All the systems above (continuous and discrete time ones) may be considered either 
on an infinite lattice (fe e Z), or on a finite one (1 < k < N). In the last case one of the 
two types of boundary conditions may be imposed: open-end (xq = oo, xn+i = — oo) or 
periodic (xo = xn, xn+i = We shall be concerned only with the finite lattices here, 
consideration of the infinite ones being to a large extent similar. 

The last remark: the "list" of references at the end of this paper might look strange; 
in fact most of the references listed in pj are relevant, and we do not reproduce them here 
only in order to save place. An interested reader is advised to consult these references. 

2 Simplest flows of the relativistic Toda hierarchy 
and their bi— Hamiltonian structure 

We consider in this section two simplest flows of the relativistic Toda hierarchy. The first 
of them is: 

d k = d k (c k - c fc _i), c fc = c k (d k+ i + c k+1 - d k - c fe _i). (2.1) 
The second flow of the relativistic Toda hierarchy is: 

d\ = d k ( k 1 j , c k = c k ( ] . (2-2) 

V^fc^fc+i d k -\d k j \d k d k+ x I 
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They may be considered either under open-end boundary conditions (d^+i = c = cn = 
0), or under periodic ones (all the subscripts are taken (mod N), so that d^+i = d\, 
c = c N , C N+ i = Ci). 

The both lattices (|1.1|) , ( |1.3|) arise from the flow ( 2A ) under two different parametriza- 
tions of the variables (c,d) by the canonically conjugated variables (x,p). Analogously, 
the both lattices (|1.2| ), (|1.3|) may be considered as arising on the same way from the flow 

(ED- 

We discuss now a Hamiltonian structure of the both flows (|2.1| ), (|2.2| ). It is easy to 
see that they are Hamiltonian with respect to two different compatible Poisson brackets. 
The first bracket is linear: 

{ck,dk+i}i = — Cfc, {cfc,c4}i = c k , {dk,d k+1 }i = c k , (2.3) 

(only the non-vanishing brackets are written down), and Hamiltonian functions generating 
the flows fl2.ip, (|2.2|) in this bracket are equal to 



i N N N 

H + = o E(4 + c fc-i) 2 + E(4 + c fc _i)c fc , H W = - ^log(4). (2.4) 

^ k=l k=l k=l 

The second Poisson bracket is quadratic: 

{cfc,Cfc + i} 2 = — CfcCfc+i, {cfc,<ifc+i}2 = —c k d k+ i, {ck,d k } 2 = Ckdk (2.5) 
the corresponding Hamiltonian functions being 

tff = E(4 + c fc ), ^ 2 ) = f:^±^. (2.6) 
fc=l fc=l 



We turn now to the integrable discretizations of the flows ( |2.1|) , (|2.2|) , derived in 
An integrable discretization of the flow ([2.1|) is given by the difference equations 



7 , a k+i — hd k+ i _ a k+ i + hc k+ i 

dk = dk n — ' Ck = Ck r~r — > ( 2 - 7 ) 

Qfc — ridk a-k + riCk 

where = Ofc(c, d) is defined as a unique set of functions satisfying the recurrent relation 

hCk--\ 

a k = 1 + hd k + (2.8) 

together with an asymptotic relation 

a k = l + h(d k + c k ^)+0(h 2 ). (2.9) 
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In the open-end case, due to c = 0, we obtain from (|2.8|) the following finite continued 
fractions expressions for a k '- 



ai = 1 + hdi] a 2 = 1 + hd 2 + 
ax = 1 + hd N H 



1 + hdi 

hc N _ x 



, , . . , h °N-2 

1 + nd N -i H 



1 + hd N n + 



+ 



ftdi 



1 + hdi 

In the periodic case ( |2.8| ), ( |2.9| ) uniquely define a^'s as A^-periodic infinite continued 
fractions. It can be proved that for h small enough these continued fractions converge 
and their values satisfy fl2.9|). 

An integrable discretization of the flow ( p.2|) is given by the difference equations 



7 , d k - ko k -i ~ c fc + hd k 

«fc — «fc+l j 7 , Cfc — Cfc + i — , l^-J-UJ 

«fc+i — Cfc + i + nQ k+ i 

where d k = ffc(c, d) is defined as a unique set of functions satisfying the recurrent relation 

Ck _ 

together with an asymptotic relation 



— = 4 - h-hok-i, (2.11) 



t> fc = ^ + 0(ft). (2.12) 

dk 



In the open-end case we obtain from ( 2.11]) the following finite continued fractions ex- 
pressions for Dfc: 



ci c 2 
*>i = 1 ri °2 



d 2 - h- 



d\ — h 

CN-l 



A U hCN ~ 

cIn—i — ft — 



— ft 

In the periodic case (|2.11|) , (|2.12|) uniquely define o^'s as A^-periodic infinite continued 
fractions. It can be proved that for h small enough these continued fractions converge 
and their values satisfy ( |2.12[) . 

It can be proved that the maps (|2.7|) , ( |2.10| ) are Poisson with respect to the both 



brackets fO), (£50 . 
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3 Lax representations 

Recall [p] that both continuous time systems fl2.1|) , ( p.2|) and discrete time ones 
( p.lO| ) admit Lax representations, the Lax matrices being the same for the both cases. 

Namely, the following statement holds. Introduce two N by N matrices depending on 
the phase space coordinates c k , d k and (in the periodic case) on the additional parameter 
A: 

N N 

L(c,d,X) = ^2 d k E kk + X^2 E k+1>k , (3.1) 

k=l k=l 

N N 

U(c,d,X) = Y. E kk-^ 1 T. c kE k ,k+i. (3.2) 
fc=i fc=i 

Here Ej k stands for the matrix whose only nonzero entry on the intersection of the jth row 
and the kth column is equal to 1. In the periodic case we have En+i,n = ^i,jv, En,n+i = 
Eff,\\ i n the open-end case we set A = 1, and En+i,n = E^n+i = 0. Consider also 
following two matrices: 

T+(c, d, A) = L(c, d, Xp-^c, d, A), T_(c, d, A) = ET^c, d, A)L(c, d, A). (3.3) 

Proposition 1. The flow (|2.1|) is equivalent to the following matrix differential equa- 



tions: 

which imply also 
where 



L = LB — AL, U = UB- AU, (3.4) 
f + = [T+,A), r_ = [T_,B], (3.5) 



N N 

A(c,d,X) = ^2(d k + c k -i)E kk + E k+lik , (3.6) 

k=l k=l 

N N 

B(c,d,X) = ^2{d k + c k )E kk + X^2E k+ljk . (3.7) 

k=l k=l 

Proposition 2. The map (|2.7[) is equivalent to the following matrix difference equa- 
tions: 

L = A l LB, U = \- l UB, (3.8) 

which imply also 

f+ = A^T+A, T_ = B _1 T_B, (3.9) 
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where 

N N 



A(c, d, A) = *kE kk + h\J2 E k +i, k , (3.10) 



fc=i fe=i 

AT N 



B(c, d, A) = J2 bkE kk + h\Y, E k+1 , k , (3.11) 
fe=i fc=i 

and t/ie quantities b k are defined by 

flfc+i — hd k+ i a k + /iCfc , . 

bfc = afc = a k-i —7 • l^.li) 

Note that the compatibility of the two expressions for b k in (|3.12|) is an immediate 
consequence of ( |2.8|) , and that it follows from (|3.12|) , (|2.9|) that 

b k = l + h(d k + c k ) + 0(h 2 ). (3.13) 



Proposition 3. The flow ( |2.2| ) is equivalent to the following matrix differential equa- 
tions: 

L = LD-CL, U = UD- CU, (3.14) 



which imply also 
where 



f + = [T + ,C], 7L = [T_,£>], (3.15) 



N Cl. 

C(c,d) = —X ^2 —^—E k)k+ i, (3.16) 
k=i 

D(c,d) = -A-^f^t+i. (3.17) 

fc=i 

Proposition 4. T/ie map (|2.10 ) is equivalent to the following matrix difference equa- 
tions: 

L = CLT>-\ U = CUF>- 1 (3.18) 

which imply also 

f+ = CT+CT 1 , f_ = DT ' D - 1 , (3.19) 

where 

N N 

C(c, d, A) = J2 E kk + hX- 1 cfc^fc,*+i, (3.20) 
fc=i it=i 

AT AT 

D(c, d, X) = E kk + hX- 1 *kEk,k+i> ( 3 - 21 ) 

k=l k=l 
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and the quantities c k are defined by 

dk - hx>k-t c k + hd k fn on v 

Ck = dk l u — = r~r • (6.22) 

— i>^k Ck+l + Mk+l 

The compatibility of the two expressions for in ( |3.22| ) is an immediate consequence 
of ( ^TTTD , and it follows from flO^ ), ( gig ) that 

c fc = ^ + 0(/i). (3.23) 



fc+i 



The spectral invariants of the matrices T±(c, d, A) serve as integrals of motion for the 
flows (|2.1|), (|2.2|), as well as for the maps fl2.7p, ( |2.10p . In particular, it is easy to see that 
the Hamiltonian functions (|2.4|), (|2.6|) are spectral invariants of the Lax matrices: 



H? = hr(Tl), F« = -trlog(T±), 

F| 2) = tr(T±), H { y = tr^ 1 ). 

Moreover, it can be proved |]J that the maps ( p.7| ), ( p,10| ) admit in both Poisson brackets 
( |2. 3D , ( |2.5|) interpolation by Hamiltonian flows with the Hamiltonian fuctions being certain 
spectral invariants of the matrices T±. 



4 Parametrization of the linear bracket 
by canonically conjugated variables 

In what follows we shall consider different Poisson maps from the standard symplectic 
space R 2N (x,p) into the Poisson space R 2N (c, d), the latter being equipped with different 
Poisson brackets, while the former always with the canonical one: 

{x k ,Xj} = {pk,Pj} = 0, {Pk,Xj} = faj- (4.1) 

We shall call such maps parametrizations of the corresponding Poisson bracket on R 2N (c, d) 
through canonically conjugated variables (x,p). 

For example, the linear Poisson bracket ( p.3|) may be parametrized by the canonically 
conjugated variables (x,p) according to the formulas: 

4 = Pk - exp(x fc - x k -i), c k = exp(x fe+ i - x k ). (4.2) 

Let us see how do the equations of motion look in this parametrization. We start with 
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Obviously, the function H+ takes the form 

l N N 

r(l) 



H+ = - J2pI + J2Pk ex P( x k+i - x k)- (4.3) 
1 k=l k=l 



Correspondingly, the flow (|2.1|) takes the form of canonical equations of motion: 
x k = dH^/dp k =p k + exp(x k+1 - x k ), 

p k = -dH^/dxk = p k exp(x k+ i - x k ) -p k -iexp(x k -x k -i). 

As an immediate consequence of these equations one gets the Newtonian equations of 
motion ( |1 . 1| ) . A standard procedure allows to find a Lagrangian formulation of these 
equations. Indeed, one has to express 

N 

C = J2ikPk~H (4.4) 
fc=i 

in terms of (x k ,x k ), which in the present case leads to 

N 

2 



1 N 2 

£+ (x,x) = - fife - exp(x fc+ i - x k )) ■ (4.5) 
1 k=i 



Note that the results of the previous section provide us with a Lax representation of our 
new lattice ( |1 . 1| ) : one needs only to set in the formulas of the Proposition 1 



c fc = exp(x fe+ i -x k ), d k = x k - exp(x fe - x k -i) - exp(x fc+ i - x k ). 
We turn now to a less straightforward case of discrete equations of motion. 



Theorem 1. In the parametrization ( f4.2|) the equations of motion (|2.7| ) may be pre- 
sented in the form of the following two equations: 

hp k = exp(x k - x k ) - r T r + hexp(x k - x k -i) -hexp(x k +i - x k ), (4.6) 

1 - hexp(x k - x k -i) 

1 

hp k = exp(x k - x k ) - 7 ztt-, (4.7) 

1 - hexp(x k+ i - x k ) 



which imply also the Newtonian equations of motion (|1.4|) . 

Proof. The second equation of motion in ( |2.7| ) together with the parametrization 
Cfc = exp(x k+ i — x k ) implies that the following quantity is constant, i.e. does not depend 
on k: 

exp(—x k + x k )(a k + hc k ) = const. 
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Choosing this constant to be equal to 1, we get: 

a k + hc k = exp(x k - x k ), 

hence 

a k = exp(x k - x k ) - hexp(x k+ i - x k ) = exp(x k - x k )(l - hexp(x k+ i - x k/ 
Substituting the last two formulas into fl2.8|), we get 

hck-i 1 



- hd k = 1 + 



or 



hd k = exp(x k - x k ) (l-h exp(x fc+1 - x k 
Now the first equation of motion in 

d k = dk 



a k -i l-h exp(x k - x k _ x ) ' 

1 



1 - hexp(x k - x k ^i) 
may be rewritten with the help of ( [4.10| ) as 

1 - hexp(x k - x k -i) 



(4.8) 
(4.9) 

(4.10) 
(4.11) 



1 - hexp(x k+1 -x k )' 



which together with (|4.22| ) implies: 

hd k = exp(x fc - x k ) (l-h exp(x k - x k -i 



1 



(4.12) 



1 - hexp(x k+ i -x k )' 

Under the parametrization dk = Pk — exp(x fc — x k _i) the equations ( |4.11|) , ([4.12|) are 
equivalent to (|^), (|77[). ■ 

The Lax representations for the system ( |1.4| ) is given by the Proposition 2, where the 
expressions for the coefficients c k , d k , ak, bk in terms of the variables Xk and their discrete 
time updates Xk are given by Ck = exp(xk+i — Xk), ( |4.11|) , (^4.12|) , ( |4.9| ), and 

b k = exp(x k - x k )(l - hexp(x k - x k - X ) 

(the last formula following from (P|), and (PHQ ). 



Note also that the equations (|4.6| ), (|4.7| ) not only immediately imply ( |1.4|) from the 
introduction, but, moreover, allow to find a Lagrangian interpretation of this equation. 
Indeed, the general theory says that if the equations of motion are represented in the 
Lagrange form 

d(A(x,x) + A(x,x))/dx k = 0, (4.13) 



then the momenta p k canonically conjugated to x k are given by 

p k = -dA(x,x)/dx k , 



(4.14) 
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so that 



p k = dA(x, x)/dx k . (4-15) 
with ( |4.14j ), (|4.15|) , respectively, we see that the Lagrange function 



Identifying 

for the equation (|1.4j ) can be chosen in the form 



N 



N 



N 



A^\x,x) = ^Lp(x k -x k )-h 1 J2^°S^--hexp(x k+ i-x k )^-J2^ x P{xk+i-x k ), (4.16) 

k=l k=l k=l 

where 

^(0 = ^ 1 (exp(0-l-e). 

Obviously, this function serves a finite difference approximation to ( |4.5| ). 

We turn now to the equations of motion (|2.2| ), ( |2.10|) , and find out how do they look 
in the parametrization ( |4.2[ ). 

The function H_ takes the form 



Hr' = - ^ log (p k - exp(x fc - 



N 



k=l 



(4.17) 



Correspondingly, the canonical equations of motion for the flow (|2.2|) takes the form: 
x k = dH®/dp k = 1 
p k = -dHr/dx h 



p k - exp{x k - x k - X , 
exp(x fc+ i - x k ) 



exp(x fc - x fc _i) 



p k+1 - exp{x k+l - x k ) p k - exp(x fc - x k ^ x ) 
As a consequence of these equations, one gets: 

1 



p k = exp{x k -a; fc _i 



Xk 



p k = -a; fc+ iexp(a; fc+ i - x k ) + x k exp{x k - 



and the Newtonian equations of motion (|1.2Q follow. A standard procedure (|4.4j) leads to 
a Lagrangian formulation of these equations. One has: 



C-Hx, x) 



N N 

log(xfc) + x k exp(x k - x k -i). 

k=l k=l 



(4.18) 



In order to get a Lax representation of the lattice ( |1.2| ) one needs only to set in the 
formulas of the Proposition 3: 



c fc = exp(x k+1 -x k ), d k = - 



1 

x k 
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Turn to the discrete equations of motion ( [2.10|) , we get: 

Theorem 2. In the parametrization ( |4.2| ) the equations of motion ( |2.1(J| ) may be 

presented in the form of the following two equations: 



h 



Pk 



exp(£ fc - x k ) 



+ exp(x k - 



V, 



(4.19) 



Pk 



exp(x fc - x k ) -1 



+ exp(x fc+ i - x k ) - exp(x fc+ i - x k ) + exp(x fc - x k -i), (4.20) 



which imply also the Newtonian equations of motion ( |1.4|) . 

Proof. The second equation of motion in (|2.10 ), rewritten as 



c k = c k 



1 + hd k /c k 

1 + /iOfc+l/Cfc+l 



together with the parametrization c k = exp^+i — x k ) implies that the following quantity 
is constant, i.e. does not depend on k: 



exp(x k -x k ) (l + — 

Choosing this constant to be equal to 1, we get: 

Cfc , , h 
h h 



const. 



exp(x fc - x k/ 



hence 



hx> k = exp(x fc+ i - x k ) - exp(x fc+ i - x k ) = - exp(x k+ i - x fe )(exp( 
The recurrent relation (|2.11| ) implies: 



d k - = — + h, 



which together with (|4.22|) , ( |4.21|) implies: 

h 



di 



exp^fc - x k ) 



- exp{x k - 



)(expf 



X k -\ — 



(4.21) 



(4.22) 



(4.23) 



(4.24) 



Now we can rewrite the first equation of motion in fl2.10|) , taking into account ( |4.23|) . 

~ exp(x fc+ i - x k+ i) - 1 
d-k — "fc+i — ; : — • 



exp{x k - x k ) 
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The last equation together with (|4.24j ) implies: 

h 



di 



exp{x k - x k ) 



1 



exp{x k+1 -x k )[ exp{x k+1 - x k+1 



(4.25) 



Under the parametrization d k = p k — exp(x k — x k -x) the equations ( |4.24| ), ( |4.25| ) are 
equivalent to ( |lg ), flOPp . ■ 

The Lax representations for the system (|1.5| ) is given by the Proposition 4, where the 
expressions for the coefficients c k , d k , D k , c k in terms of the variables x k and their discrete 
time updates x k are given by c k = exp(x fc+1 — x k ), (|4.24|) , (|4.25|) , (|4.22| ), and 



h 



c k = -exp{x k+1 - x k )\exp{x k+ i - x k+i 



(the last formula following from flOg ), Q), and (|421|) ). 

Identifying ( |4.19| ), ( |4.2(J| ) with ( [4.14[ ), ( ^4.15| ), respectively, we get the Lagrange function 
for the equation ( |1.5| ) in the form 



N 



N 



X, X 



-h ^( x k -%k) + J2 [ ex P( x k ~ x k -x) ~ exp(x fc - x k _x) 
k=i 



k=l 



where 



^(0 



drj 



log(exp(f) 



io exp(r]) — 1 

This function clearly is a finite difference approximation of ( [4.18| ). 



5 Parametrization of the quadratic bracket 
by canonically conjugated variables 

We give for completeness the results corresponding to another parametrization of the vari- 
ables Cfc, d k by means of canonically conjugated variables x k ,p k , namely the parametriza- 
tion leading to the quadratic bracket ( |2.5| ). The relativistic Toda lattice arises on this 
way. The corresponding formulas were given in jl], but in an ad hoc manner, without 
derivation. We take an opportunity to fill in this gap here. 

The parametrization leding to the quadratic bracket (|2.5|) reads: 

4 = exp(p fc ), c k = g 2 exp(x fc+ i - x k + p k ), (5.1) 



(g 2 G 1 is a coupling constant) 
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In terms of these variables 

N 



#j 2) = ^exp(p fc )(l +g 2 exp(x k+1 -x fc )), 

k=l 

N 

H {2) = exp(-p fe )(l + g 2 exp(x k - x k -i)). (5.2) 



k=l 



Hence the equation of motion corresponding to ( |2.1| ) take the canonical form 

x k = dH^/dp k = exp(p k ) (l + g 2 exp(a; fc +i - x k )) , 

p k = -dH { + ] /dx k = g 2 exp(x k+1 - x k + p k ) - g 2 exp(x k - x fe _i +p k -i). 

This can be put into a Newtonian form (|1.3|) . 

A standard procedure ( [4.4| ) allows also to find a Lagrangian formulation of these 
equations, with a Lagrange function 

N N 

£ { +\x,x) = J2[x k log{x k ) -x k ] - ^x k \og (l +g 2 exp(x k+1 -a; fc )). (5.3) 

k=l k=l 

The Lax representation for these equations are given by the Proposition 1 with 

dk = %k/ (l + 9 2 exp(x fe+ i - x fe )) , c fe = g 2 exp(x fc+ i - x k )d k . 

Analogously, the canonical equations of motion corresponding to ( |2.2| ) are: 

x k = dHr/dpk = -exp(-pk)(l + g 2 exp(x k -x k -ij), 

p k = -dH^/dxk = g 2 exp(x k+ i - x k - p k+ i) - g 2 exp(x k - x fc _i -p k ). 



The Newtonian equations following from these ones are just the same ( |1.3| ) as before, 
they correspond, however, to a different form of a Lagrange function: 



N N 



Cl)(x,x) = - ^[x fe log(-x fe ) - x k ] + ^fclog (l + g 2 exp(x k - x k -i)). (5.4) 

k=l k=l 

Respectively, the Lax representation for these equations is given by the Proposition 3 
with the identifications 

l + g 2 exp(x k -x k -i) 2 , 

dk = : , c k = g exp(x k+1 - x k )d k . 

Xk 



We turn now to the discrete time systems (|2.7p , ( |2.10| ). 
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Theorem 3. In the parametrization ( 573]) the map (|2.7| ) tofces the form of the following 
two equations: 

(exp(x k -x k ) - l) (l + # 2 exp(x fc - x k -i)) 

hexp(p k ) = -j (5.5) 

1 + g 2 exp(x k+1 - x k )) (1 + g 2 exp(x k - x k -i) 

exp(x k - x k ) - 1 



hexp(p k ) = jf ^-r- (5.6) 

(^1 + ^exp(x fc+ i - x fc ) J 

T/ms implies also a Newtonian form (|1.6|) o/ equations of motion. 



Proof. From ( |2.7p it follows: 

Cfc _ c k a k+ i + hc k+ i a k — hd k 
d k d k a k+ i — hd k+ i a k + hc k 

Since c k /d k = g 2 exp(x k+ i — x k ), this implies that the following quantity is constant, i.e. 
does not depend on k: 

,~ \ a k — hd k 

exp(xfe — Xk) : — = const. 

a k + hc k 

Setting this constant equal to 1, we get: 

a k + hc k 



This implies: 



a k 



. , exp(x k - x k ). 
a k — na k 

d k exp(x k - x k ) + c k 
exp(x k -x k ) - 1 

exp(x k -x k )(l+g 2 exp(x k+1 -x k ) 

hd k ^ r '-. 5.7 

exp(x fc -x k )-\ 



h 



As a consequence, we get: 

d k + c k 



a k + hc k = hexp(x k - x k ) 



exp(x k - x k ) - 1 



exp(x k - x k ) (l + g 2 exp(x fc+ i - 



= hd k K - r— : '-■ (5.8) 

exp(x fe - x k ) - 1 

Substituting (p.7|), ( |5.8|) in the recurrent relation ( |2.8| ), we get: 

, , - , hc k _x 1 +g 2 exp(x k - x k -i) . . 

ak ~ hd k = H = — — = -, — r. (5.9) 

a fc _i 1 +g 2 exp(x k - x k -i) 
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Substituting in the left-hand side of this formula the expression ( |5.7|) for o^, we arrive at: 

(exp(x fc -x k ) - 1) (l + g 2 exp(x k - x fc _i)J 



hd k — 

(l + g 2 exp(x k+1 - x k )J (l + g 2 exp(x k - x k -i) 
Further, from the first equation of motion in Q2.7| ) and ( |5.9| ) it follows: 

exp(x k - x k ) - 1 



(5.10) 



hdk 



1 + g 2 exp(x k+ i - x k ) 



(5.11) 



Now ( |5.5| ), ( |5.6| ) follow from ( |5.10| ), ( [5.11| ) under the parametrization d k = exp(p k ). 



The Lax representation for the system ( |5.5|) , (|5.6|) is given by the Proposition 2 with 
the following expressions through x k , x k . ( p.lO[ ), ( [5. 1 1| ) for d k , d k ; c k = g 2 exp(x k+ i—x k )d k ; 



<*fc 
i>k 



exp(x fe - x k 



(l + g 2 exp( 



x k+l ~ x k 



)) (l + g 2 exp(x k - 



X k -1 



(l + g 2 exp(x k+ i - x k j"j (l + g 2 exp( 



x k - X k -i 



exp{x k - x k ) 



(5.12) 
(5.13) 



Indeed, ( |5J3D follows from ( ^TUQ , and (^Wj ) follows from QgTTg ), and 

Identifying (|5.5| ) and (|5.6|) with ( |4.14| ), ( |4.15|) , we get a Lagrange function 



N N 

A$\x, x) = J2 -%k) + i^( x k+i - %k) - ^(x k+ i - x k )} 

k=l k=l 



(5.14) 



where the two functions $(£), ^(0 are defined by 



log 



exp(r/) — 1 



*(0 = J o log (l + g 2 exp(r/)) c?77- 



(5.15) 



It is easy to see that this Lagrangian function serves as a finite difference approximation 
to© 

Theorem 4. In t/ie parametrization (|5.1|) i/ie map (|2.10| ) iafces £/ie /brm o/ t/ie follow- 
ing two equations: 



exp(p k ) 
exp(p fc ) 



h(l +g 2 exp(x k -x k -i 



(l - exp(x k - x k 
h(l + g 2 exp(x k - x fc _i 

(l - exp(x k - x k )) (l + g 2 exp(x k+1 - 



1 +5f 2 exp(x fc+ i -x k \ 



x k . 



(5.16) 
(5.17) 
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This implies the same Newtonian equations of motion (|1.6|) as for the map (|2.7f) . 
Proof. From ( fUCQ , (g7TT|) we deduce: 

dfc ^fc+i dk — hdk-i c k+ i + hOk+i dk+i <>fc+i 

Because of Ck/dk = g 2 exp(xk+i — Xk), this implies that the following quantity does not 
depend on k: 

d k exp(x k - Xk+i) = const. 
Setting this constant equal to g 2 , we get: 

,2 



9 fc = g exp(x k+ i - x k ). 
Substituting this formula into the recurrence ( |2.11|) , we get: 



hll + g 2 exp(x k - x k -\ 



Formulas (|5.16|) , (|5.19Q imply also 



dk — hd 



1 - exp(x fc - x k/ 



h(l + g 2 exp(x k - x k -i 



(5.18) 



(5.19) 



fc-i 



1 - exp(x fc - x k ) 

This last formula together with the first equation in ( |2.10|) implies: 

h(l + g 2 exp(x k - x fc -i)) (l + g 2 exp(x fc +i -x k ] 



(5.20) 



d k 



(l — expl 



00 fa 00 fa 



1 +5f 2 exp(x fc+ i - x k ] 



(5.21) 



Finally, ( |5.16| ), ( |5.17| ) is equivalent to ( |5.19| ), ( |5.21| ), due to the parametrization d k = 
exp(p fc ). ■ 

The Lax representation for the system ( (5.16j ), (|5.17|) is given by the Proposition 4 with 
the following expressions through x k , Xk'- fl5-19|) , (|5.21|) for d k , d k ] Ck = g 2 exp(x fc+1 — x fc )(4; 
Q5.18Q for d k ; and 



c k = g exp(x k+ i - x k 



(l — expl 



Xk+l ~ x k+l 



1 — exp(xfc — Xk] 



1 + g 2 exp(x k - x k -\] 
1 +g 2 exp(x k +i - x k ] 



The last formula follows from ( ggg ), (^T§, and ( ^U|) . 

Identifying (|5.16|) and ( |5.17| ) with ( |4.14j ), ( |4.15| ), we get a Lagrange function 



A {2) (X,X) = -J2 ®( X k ~ X k) + Yl ~ X k-l) ~ ^( X k ~ x k-l) 



N 



N 



(5.22) 



k=l k=l 

wit the same functions $(£), (|5 . 1 5 ) as before. It is easy to see that this Lagrangian 
function is a finite difference approximation to ( |5.4|) . 
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6 Parametrization of the mixed brackets 
by canonically conjugated variables 

It turns out that there exist still another parametrizations of the variables (c, d) by canon- 
ically conjugated variables (x,p) leading to interesting discretizations. As we shall see, 
these parametrizations lead to Poisson brackets which are linear combinations of the two 
homogeneous ones fl2.3p and ( |2.5| ). In some sense (which will be clear from the proofs of the 
Theorems below) these parametrizations are specially designed to obtain nice Newtonian 
equations from the maps (|2.71), ( j2.10| ). 

We start from the parametrization leading the linear combination {•, -}i + h{-,-}2, 



which turns out to admit nice discrete Newtonian formulation when applied to 
Clearly, we get an alternative discretization of the lattice on this way Consider 

following parametrization: 

hdk = exp(hp k ) - 1 - hexp(x k c k = exp(x k+1 - x k + hp k ). (6.1) 



(Obviously, in the limit li^Owe recover the parametrization of the linear bracket (|4.2j) ). 
Simple calculations show that the Poisson brackets between (c, d) induced by ( |6.1| ) read: 

{c k +i,c k } = hc k+ ic k , {d k+ i,d k } = —c k , 
{d k +i, c k } = c k + hd k+ ic k , {d k , c k } = —c k — hd k c k , 

which is exactly {-, -}i + h{-, -}2- Let us look at the equations of motion generated by 
these /i-dependent Poisson bracket. 

Theorem 5. In the parametrization (|6.1|) the map (|2.71) takes the form of the following 



two equations: 

exp(hp k ) = exp(x k - x k ) ± ^—j- v ^r— ' L ^ J - , (6.2) 

(1 + hexp(x k+1 - x k )J 

exp(hp k ) = exp(x fc - x k )(l - hexp(x k+1 - x k j). (6.3) 



(l + hexp(x k - x k -i)J 


(l - hexp(x k - x k -i)j 


(l + hexp(x k+1 - x k )^ 





This implies the Newtonian equations of motion (|1 . 7|) . 

Proof. The crucial observation consists in the following: we can extract from the 
relations ( |6Tl| ) the following consequence: 

exp(hp k ) = 1 + hdk H , ■ 

exp(hpk-i) 
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Comparing this with (|2.8| ), we see that and exp(hp k ) satisfy one and the same recurrent 
relation. Due to the uniqueness of its solution, they must coincide, so that we obtain: 

a k = exp(hp k ). (6.4) 

As a consequence, we get immediately: 

a fc - hd k = 1 + hexp(x k - x k -i), a k + hc k = exp(hp k )(l + hexp(x k+1 - x k )J. (6.5) 

These expressions together with ( |6.1| ), being substituted into (|2.7|) , allow to rewrite the 
latter in the form: 

exp(hp k ) - hexp(x k - x k -i) = exp(hp k ) — — - h exp(x k+1 - x k ); 

1 + hexp(x k - x k -i) 

ew(lM _ ~ xt + m = eMxt+1 _ xt + hpt+l) l + ftexpfa. +2 -. i+1 ) (66) 

1 + hexp{x k+ i — x k ) 

The formula arising when excluding exp(p k ) from these two equations, can be we written 
after some manipulations as: 

/ — , 7 l + hexp(x k+2 -x k+1 ) _ 

exp(a; fc+ i - x fe+ i + hp k+1 ) — + hexp(x k+1 - x k ) = 

1 + n exp{x k+ i — x k J 

i ~ . ,l + hexp(x k+1 -x k ) „ 

exp(x fc - x k +p k ) — — r + hexp(x k - x k -i). 

1 + hexp(x k - Xk-i) 

So, the expression on the right-hand side is constant, i.e. does not depend on k. Setting 
this constant equal to 1, we arrive at (|6.2|). Substituting ( |6.2|) into ( |6.6| ), we get ( p.3| ). ■ 
It is not difficult to extract from this proof the expressions for coefficients of the 
matrices forming the Lax representation of the system ( |1.7| ) following from the Proposition 
2. Also a Lagrangian formulation of this system can be obtained in a standard way: the 
Lagrange function corresponding to fl6.2p, fl6.3| ) is 

d N fjjj k _ Xk \2 N 

A^ llxcd) (x,x) = k n h k h l 12 ~ x k ) + 4> 2 {x k+ i - X k ) , (6.7) 

k=l Zl1 k=l 

where 

<M0 = J lo S ( x - hexjp(r]))dr}, 2 (O = J log (l + hexp^dij. 
This is a finite difference approximation to (|4.5| ) different from (|4.16|) . 
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We turn now to another parametrization leading to a mixed Poisson bracket, namely 
to the bracket {•, — h{-, -}\. The corresponding formulas are: 

d k = exp(p k ) + h(l + g 2 exp(x k - x k -ij), c k = g 2 exp(x fc+ i - x k +Pk)- (6.8) 
As is easy to calculate, the resulting Poisson brackets between the variables (c, d) are: 

{cfc+i,Cfc} = Cfc+iQc, {dk+i,dk} = —hck, 

{4+i, c fc } = 4+iCfc - hc k , {d k , c k } = -d k c k + hc k . 

i.e. the linear combination {•, -} 2 — h{-, - }\. The equations arising from (|2.10|) under this 
parametrization, naturally approximate the relativistic Toda lattice ( |1.3j ). 

Theorem 6. In the parametrization (|6.8|) the map (|2.10| ) takes the form of the follow- 
ing two equations: 



h(l + g 2 exp(x k — x k _i) 

exp{p k ) = —, — , (6.9) 

yexp{-x k + x k ) - 1) 

h(l + g 2 exp(x k - x k -i)) (l + g 2 exp(x k+1 - x k ) 

eMPk) = 7 ( ~ I ^ A h I 2 F ^ ((U()) 

(^exp(-x fc + x fc ) - 1J (^1 +5f 2 exp(x fc+ i 



T/izs implies the Newtonian equations of motion (|1.8|) . 

Proof. This time the crucial observation consists in the following: as a consequence 
of the relations fl6.8|) we have: 

— r = d k - h - hg 2 exp(x fc - x k -i). 

g 2 exp(x k+1 - x k ) 

Comparing this with fl2.11|) , we see that d k and g 2 exp(x k+ \ — x k ) satisfy one and the same 
recurrent relation. Due to the uniqueness of its solution, they must coincide, so that we 
obtain: 

D k = g 2 exp(x k+1 - x k ). (6.11) 
As a consequence, we get immediately: 

d k - ho k -\ = exp(p k ) + h, c k + hx> k = g 2 exp(x k+1 - x k )(exp(p k ) + fij. (6.12) 

Substituting these expressions together with ( |6.8| ) into ( |2.10| ), we bring the latter to the 
form: 

/ — \ I i 2 /- ~ v / \ I * 2 / \ exp(p fc ) + h 

exp(p k ) + hg exp(x k - x k ^) = exp(p k ) + hg exp(x fe+ i - x k ) 



exp{p k+1 ) + h 



,~ ~ . ~n / . \ exp p fc + h 

exp(x k+1 - x k + p k ) = exp(x k+ i- x k +p k+1 ) r— — . (6.13) 

exp(p k+1 ) + h 
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Excluding exp(p k ) from these two equations, we get the formula which after some manip- 
ulations can be written as: 

exp(x fc+ i - x k+1 +p k +i) - hg 2 exp(x k+ i - x k ) _ exp(x k - x k + p k ) - hg 2 exp(x k - x kl ) 
exp(p k+1 ) + h exp(p k ) + h 

So, the expression on the right-hand side is constant, i.e. does not depend on k. Setting 
this constant equal to 1, we arrive at (|6.9D, which is equivalent also to 



hexp(-x k + x k )(l + g 2 exp(x k -x k „i] 

exp(p k ) + h = v '-. 6.14 

exp(-x k + x k ) - 1 

Finally, substituting ( |6.14|) into ( |6.13| ), we get 

exp(-S k+1 + x k+1 )-l) (l+fewfrSk-i 



exp(p k ) = exp(p t 



k+l, 



(exp(-x fc + x k ) - l) (l + g 2 exp(x k+1 - x k ] 



which together with (|6.9|) implies ( |6.10| ). I 

The Lax representation for ( |1.8j ) is given by the Proposition 4; it is not difficult to 
extract from the proof above the expressions for the entries of the matrices forming the 
Lax representation. Also the Lagrangian formulation of (|1.8[) easily follows from ( |6.9|) . 
( |6.10| ). The corresponding Lagrange function is 



N N 



A (mixcd \x, x) = j2 $ (-^ + **) + E - 5 *-0 - - ^-0] » ( 6 - 15 ) 



fc=l fc=l 



wit the functions ^(0 given in (|5.15 ). 

Let us note that, though the equations flOD , (|1.8|) are very similar, there exist no 
obvious changes of variables bringing one of them into another. The only way to do this, 
i.e. to connect two corresponding sets of (x k , x k ), is to identify the corresponding (c k , d k ), 
given for flL|) by (§!]), ( ^T6|) , QOTQ , and for (|T|) by (|B^), ( |6TT0|) . The resulting 

change of variables is a rather nontrivial Backlund transformation. 



7 Conclusion 

The main message of the present note is following: the field of integrable systems of 
classical mechanics, even in its best studied parts, is far from being exhausted. Namely, 
the well known flows of the relativistic Toda hierarchy fl2.1|) , (|2.2j ) have a much more 
rich dynamical content as usually assumed. Even more is this true for the recently de- 
rived discretizations (|2.7|) , ( |2.10|) of these flows. Namely, different parametrizations of 
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the variables (c,d) by canonically conjugated variables (x,p) (corresponding to the bi- 
Hamiltonian structure of the relativistic Toda hierarchy) allowed us to derive two new 
integrable continuous time lattices and four new integrable discretizations, in addition to 
the previously known ones. All the obtained systems are connected by means of highly 
nontrivial Backlund transformations, which can be obtained by identifying the variables 
(cfc, dk) for different models corresponding to one and the same (continuous or discrete 
time) flow flOp, O, flT7D, or CT) . 

The methods of this paper can be used also in a more simple situation of the usual 
Toda lattice, where they also lead to interesting findings. They will be reported in a 
separate paper. 
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